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I. SUMMARY 


The logic of the subroutines used by the Solar and Lunar Orbit Perturbation 
Effects computer program, SLOPE, along with their relationships to the main 
program and to each other, is presented. 
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H. INTRODUCTION AND EXPLANATION OF THE LOGIC CHARTS 


The Solar and Lunar Orbit Perturbation Effects computer program, 
SLOPE, is designed so that the basic computations are performed by several 
subroutines. This paper presents the logic of these subroutines and describes 
their relationships to each other and to the main program. It is intended to be 
used as an aid in correlating the computations to Reference 2 and to show the 
logical flow of these computations. 

The Common statements, Data statements and Format statements for 
each subroutine are presented along with a Table of Contents and Cross- 
Reference Listing for the subroutines. The Cross-Reference Listing show 
for each label in the subroutine the page and symbol number for the label 
itself and for all points in the flow charts at which transfers to it are found. 
These page numbers refer to the flow charts and are at the right uppermost 
comer of each flow chart. They are not to be confused with the page numbers 
of this paper itself. An asterisk, *, indicates that there are other transfers 
to a label other than the one shown. 
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HI. MAIN PROGRAM, SATD 


The main program does a minimum of computation. Its main function is to 
control the several computational options and to control the input and output of 
the program. The logic charts presented on the following pages correlate the 
interaction of the subroutines and the main program. (See page 2 of Reference 
1 .) 



COLE SATD 


DECK $AT0 

COMMON/ORB /Q (28) ,A,E,0l,EL,W,CW,EN,I6UG,L,M 
COMMON /DER/D (6) , WO , C WD , E LD , 1 SE C ( 2 ) ,KSL 
COMMON/RAO/P (6,6) ,B(6,6) 

COMMON/PER T/ARG (125 ,2) ,ARGD (125,2) , PAR (125,2,6) ,SEC(5,2) 
COMMON/COEF/C (125,5) 

COMMON/CONS/ C J 2 , C J 4 , AE , PL I M 

DIMENSION SAT(12) , A SL ( 3 , 3) ,1C(5) , KSL ( 2) ,PT(5) , PG ( 5 , 2) , TOP ( 5 ) 
DIMENSION STA (5) , TERM (5) 

DIMENSION 01 OUT (6,10) , I C 1 (6,3) ,PARD1 (6,6) , FLG 1 (6) 

DIMENSION DUMOUT (123,13) , ICO (123, 5) , P ARDUM ( 1 23 , 6) ,FLGDUM(123) , 

I TRMD (123) 

DIMENSION BLOCM1QO) 

DIMENSION MESS(4,5) , A $ ( 4 ) 

EQUIVALENCE (DUMOUT , D10UT) 

EQUIVALENCE (D10UT,IC1) , (D10UT(1,4) ,PARD1) , (D 1 OUT ( 1 , 10) , FLG1 ) 

V 

EQUIVALENCE (DUMOUT (1 , 1 ) , 1 CD ( 1 , 1 > ) , (P ARDUM ( 1 , 1 ) , DUMOUT (1 , 6) ) , 
(FLGDUM(l) , DUMOUT (1,1 3)) , ( I TRMD ( 1 ) , DUMOUT (1 , 12) ) 

EQUIVALENCE ( I SEC ( 1 ) , I SI ) , (I SEC (2) , I $2) 

DATA (ASL (I , 1) , I = 1,3)/ 6H SOLAR , 6H&R A V I T , 5H A T I ON / , ( A SL ( I , 2) , I = 1 , 3 ) 
/6HLUNAR ,6HGRAVIT,5HATI0N/ , ( A SL ( I , 3 ) , I = 1 , 3) / 6HS0L AR ,6HPRE$SU, 
2HRE / 

, TWP/6. 2831853/ , VUPD/ 107. 085/ , DPR/ 5 7. 29 578/ ,RPD/ .0174532925/ 

, FLAG / .99999999EOO/ ,0NE/1 .0/ 

OATA (MESS(1,I) ,I = 1,5)/6NPE 0RB,6HIT GEN , 6HER A TOR , 6H WAS U,3HSED/, 
(MESS(2,I) ,1=1,5) /6HGILL 0,6HRBIT G , 6HENER AT , 6HOR WAS,6H USED /, 
(MESS(3,I) ,1=1 ,5) /6HBRWER ,6H0RBIT , 6HGENER A , 6HT0R WA,6HS USED/ , 
(MESS (4,1) ,I=1,5)/6HHST 0R,6HBIT &E , 6HNERAT0, 6HR WAS ,6HUSED / 

DATA ACUT/3HCUT/, A0LNK/3H / 

DATA (ASCI) ,1=1,4) /6H ,6HS0LAR , 6H , 6HLUNAR / 

EQUIVALENCE (K$ , KSL ( 1 ) ) , (KL , KSL (2) > 

1000 FORMAT (12A6) 


1010 FORMAT (4E12. 8) 

1020 FORMAT (4 1 6, F6 • 0) 

1030 FORMAT (3F12*6,2I6,2I3,316) 

1040 FORMAT (F6.0, IX, F2.0, I 2 , F2 . 0 , 1 X , 2F2 . 0 , 1 X , F2 . 0 , F3 . 0 ) 
1050 FORMAT (6E12 • 8) 

1100 FORMAT (F10.5, 1X711) 

1500 FORMA T ( 5E 14 • f ) 

2000 FORMAT ( 1 Ml // 30X2THMM* START OF CASE 6**8*) 
5010 FORMAT ( 1H0//3X1 2A8/) 
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COll I* TO 


(COMT) 


MISCELLANEOUS statements 


OIC A SATO 

SOU FORMAT ( 1 H 1 , S X i 1 2 A 9 ) 

2010 FORMAT ( 1M0SX1 3HEPOCH DATE ISF11.8,10H CENTUR 1 E 5 , 20X , I 2 , 1H / , 1 2 , 1 M / , 
12.F12.6 // 6X , 2 1 H SA TELL 1 TE HAS AREA OF E14.8,Z0H SO. CM. AND MASS 
OF Ell .4 ,8H CRAMS ) 

2020 FORMAT UM010X3A6) 

2050 FORMAT (1H010X5H0AY F7.2.10H OF EPOCH /// 

7X1THPERTURBINC FORCES) 

2040 FORMAT (1H010X29HSOLAR RADIATION PRESSURE, F =E15.B,14H CM./SEC./SE 
C.) 

2050 FORMAT (///20X68HORBITAL ELEMENTS (EXCLUOING LUNAR-SOLAR PER TURBA T I 
ONS AT THIS POINT) 

//12X1HA19X1HE19X1HI19X1HL/ 
4E20.8//9 X 5HOMEO A 1 6X9HC AP OMEC A 1 3 X 1 HN 1 9 X 1 H T / 4E 20 . 8 / / ) 

2080 FORMAT (1H0,€X28HARGUMENT$ OF TERMS HAVE FORM// 

10X81H (N14LAMBDA PRIME ♦ N2*0MEGA PRIME ♦ N34CAP OMEGA PRIME ♦ N4 
♦OMEGA ♦ N5*CAP OMEGA)//) 

2090 FORMAT (1H012X26H INCLUDING SECULAR TERMS - 2A6 ) 

2100 FORMAT (26X7HDELTA E9X7HDELTA I 9 X 7HDEL T A L7X1 1HDELTA OMEGA 
3X15HDELTA CAP OMEGA/) 

2110 FORMAT ( 2X I 4 , 2X 5 I 3 , 6E 1 6 . 8 , I 4 , 1 X A 3 ) 

2120 F OR MAT (11X313, 6E 16. 8, IX, A3) 

2200 FORM AT (1H014X5HTOTAL,5E 16.8, 3X,3A6 ) 

2210 FORMAT(1H012X7H$ECULAR,5E1€.8,3X,3A6 ) 

2300 FORMAT (1H08XHHGRAND TOTAL , 5E 1 6 . 8 ) 

2310 FORMAT (9X14HN1 N2 N3 N4 N5,6X 

7HDELTA E9X 7HDEL T A I 9 X 7HDEL T A L7X 1 1HDELT A OMEGA 
3X1 5HDELTA CAP OMEGA , 3X , 6HPER I OD , 6X4H TERM) 

2320 FORMAT (12X8HN4 N5 N1,6X 

7HDELTA E9X7HDELTA I 9X 7HDEL T A L7X 1 1HDEL T A OMEGA 
3X15HDELTA CAP OMEGA , 3X , 6HPER I OD ) 

2400 FORMAT (1H028X1HE , 15X1HI , 15X1HL , 13X5HOMEGA .9X9HCAP OMEGA// 
20X5E16.8) 

2420 FORMAT (1H0, 12X, 1 HA , 1 5 X 1 HE , 1 5X 1H I , 1 5X 1 HL , 1 3X 5H0ME G A , 9X9HC AP OMEGA/ 
/ 4X6E16.8) 

2410 FORMAT (1H0.3X, 16HPERTAPE ELEMENTS 5E16.8 ) 

2500 FORMAT (1H010X21HAMPLI TUDES BY TERM - 3A6//) 

2510 FORMAT (1H010X36H AMP LI TUDES BY TERM - SOLAR RADIATION//) 

2600 FORMATUH010X16MC0NSTANTS * J2 =E15.8 l 6H t J4 =£15.8, 

1 4H , PER. LIMIT =E 1 5 . 8 / 20X , 1 6H SOL AR PRESSURE =E15.8 ,13H DYNES/S 
Q • CM. ) 

3000 FORMAT (1H120X12HDEBUG OUTPUT/) 

3100 FORMAT(lH0 t 10X,20HSEMl-MA JOR AXIS, A =F12.8) 

3200 FORMAT(1HO,10X,1 7HECCENTR I C I TY , E =F12.8) 

3300 FORMAT ( 1H0 , 1 OX , 48HRED UND ANC Y OR END OF FILE ON BINARY ELEMENT TAPE 
) 
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COll UTD 


DICK 


S A TO 

3400 FORMAT (1H0,10X,2BHBINART ELEMENT TAPE INPUT - 3Af ) 

3500 FORMAT (IHO ,10X ,47H0EITA T OPTION MUST BE USED WITH PERTAPE OPTION) 
3700 FORNAT<1HO»12X3QH** WARNING - A OUT OP RANGE **//> 

3900 FORNATUHO.IBX^BH*** ERROR STOP ***) 

3910 FORMAT ( 76H 1 NO PERTURBATION OPTIONS SELECTED ON TIME ANO PERTUR 
BAT ION OPTION CARO NO. M » IBM - CARO IGNORED ) 

3920 FORMAT (4BH1 TIME ON TIME AND PERTURBATION OPTION CARO NO- 14 , 

47H IS BEFORE TIME ON PREVIOUS CARD - CARD IGNORED > 

3930 FORMAT ( 78H1 NO PERTURBATION OPTIONS SELECTED ON PERTURBATION OP 
TION CARD - CASE IGNORED ) 
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wR I ’C 

BCD 

UNI ’ 

M 

forma t 

2000 


(BLOCK , IN) 


AZ z YEAR * 363. i 
A 1 N T ( (YEAR-1 .) /A.) 


WRITE BCC 

UNIT M 

FORMAT 5010 


COMPUTED GO TO 
FOR IN 

15 AT 1.33 

920 AT 9.06 
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I/O LIST IS 
SAT 



1 1 


BRANCH ACCORDING TO THE 

VALUE of variable in 


WRITE BCC 
UNIT M 
FORMAT 2600 


I/O LIST IS 
CJZ, CJ4, PLIM, 
SOLP 


RE AC BCD 

uni; L 

FORMAT 1010 




14 

□ 

I/O LIST IS 

AREA f EM 
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BEGIN DO LOOP 
161=1, MONTH 
; 

SET I = 1 



2.04 

1 .32 —+ 

12 
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COMPUTED GO TO 


FOR 

I 



1 1 

AT 
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. 36 
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AT 
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AT 
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. 37 
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AT 

1 
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READ BCD 
UNIT L 
FORMAT 1020 


BRANCH ACCORDING TO THE 
value OF VARIABLE I 








































coit IATC 


OtCA SMC 


M(( t 

(CONT) 
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COLE SATD 


DECK SATD 


DELTP = OELTl 


(TENC-OELTO) /OElTl 
♦ l . 

CELT 1 = 

DEL TI / 3 65 2 5 . 

T = TZ ♦ 

DEL TO/36525 . - 

DEL T 1 


I DTO .NE. 0 


IOFW .EQ. 0 


READ BCD 
UNIT L 
FORMAT 1100 


T = T - 3. 

* DELTI 

K T = K T 

♦ 3 

PA = 

A 

PE = 

E 

POI = 

01 

PL = 

EEL 

PW = 

£W 

1 PCW = 

j 

EC W 

FK T 5 

K T 

PER = TWP/EN * 
I 3 # 44 72 



I/O LIST IS 
ET , KSL ( 1 ) , K$L (2) 
1 SEC , KR, 10, IBl 





•5 011 I* TO 


13.4472 IS 
CONVERSION factor 
FRO* CANONICAL UNIT 
or TIME TO 
wiNurcs 


CAYS = 

CAYS ♦ 1 . 

fnon = 

MON T ♦ 1 

CELT© = 

• 

CCLTO - 3. 

DELTP 

001 = 

01 * DPR 

DC W = 

ECW * C PR 

CW = 

EW * DPR 

del = 

EEL * DPR 



/ WRITE BCD 
/ UNIT M 
f FORMAT 3920 


1 

1 

I/O LIST IS 

17 

~1 

6 


J 


REWIND 
UNIT ITF 


WRITE BINARY 
UNIT I TP 


I/O LIST 

I S 

OELTP, FMON, 

DAYS, 

YEAR, SATID, 

A , E , 

DOI , DC W 

CW 

, DEL , 

TEFM , PER 

FKT , ONE | 


BEGIN DO L OOF 


OtCK IATC 


F A 4 f i 

(CONT) 


Xsi !1 ) ♦ * DL (2^ X ' ♦ ) 

< + T * I SEC f 1 ) / 

♦ I SEC 2) X 


24 19 

/ 7 

/ WRITE BCD / 
/ UNIT V / 


T - 12* ET/3652 


DE.T = ET * VUFC 
DOT 
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WRITE BCD 
UNIT M 

FORMAT 1 5 1 0 















COLI ’SATO 


DtC« SMC’ 



T 


BEGIN CO LOOP 
TO J = 1 , 6 
; 

SET J = 1 


% 


/ 

8 

p ( J , 
TWP/ (B ( J , 

, 6) = 

,4) *VUPD) 



HIGH 




P * G C 5 


(CONT) 
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WRITE BCD 
UNI T M 
FORMAT 2200 


I/O LIST IS 
<PT (K) , K*1 , 3) , 

<A$L (I J,3> ,1 J*1 ,3) 


© 


KR .ne. 0 


DL = 

AMOD ( PL ♦ TOP (3) *RP0* 
TWP, TWP > 


QW = 

AMOD tPW* TOP ( 4 ) *R P 0 ♦ 
TWP » TWP ) 


DC W = 

AMOO (PCW^TOP (3) *RPD 
TWP, TWP) 


WRITE BINARY 
UNIT l TP > 


0 

















MtC • 

(CONT) 































IV. SUBROUTINE SIX 


This subroutine computes the set of six derivatives: 

dl dl do> dw dQ d£2 

de ’ di ’ de ’ di ’ de ’ di 

from pages 8, 9 and 43 of Reference 2. (See page 22 of Reference 1.) 
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DECK SlXD 

COMMON/ORB/ 9 (28) , A , E , 01 , EL , W , C W , EN , I BUG , L , M 
COMMON/DER/D (6) ,WD,CWD,ELD 
COMMON/CONS/ C J 2 , C J 4 f AE , PL I M 

3100 FORMAT (1H010X9H0MEGA DOT , El 6 . $ , 5X 1 3HC AP OMEGA DOT , E 1 6 . 8 , 5 X 5HL DOT, 
E16.8 , 3X 1 1 2H (OBLATENESS) ) 

3200 FORMAT (1H03X11HD(L DOT) /DE.5X1 1HD (L DOT) /DI , 3X15HD (OMEGA DOT)/DE, 

1 X 1 5HD (OMEGA DOT ) /D 1 , 2X 1 *HD ( C AP OMEGA DOT ) / DE , 2X 1 9HD ( CAP OMEGA DOT 
) /DI/4E16.8 , 2E20.8) 
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<<Kt UTO 


DEC* sue 


Pktl 10 


JUMOUTINC Ml 




T 5 = - 15. * A J 4 * 

E2 * 

(5. -30.*C0UO5.* 
COM) /DEN 





\JC0NT) 

Q 


41 
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COLE SATO 


DECK SI XD 


M&E 11 


subroutine six 
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V. SUBROUTINE SUN 


Subroutine SUN computes solar ephemeris quantities and their derivatives 
with respect to time using the relationships on page 46 of Reference 2. These 
quantities are referenced to the earth’s equatorial plane. (See page 22 of 
Reference 1.) 



DECK SUNO 


COMMON/ORB/Q$<7) , DQS t 7) , QL ( 7 ) , DQL ( 7) , A f E , 01 , EL , W , CW , EN 
DIMENSION C f 1 5 ) 

DATA (C(I),I = 1,151/0. 999997,0. 01675104, - 0 . 4 1 8E-4 , -0 . 1 3E - 6 , 

0 . 98 560027 , 279 . 69 668 | 36000 .7689 , 0. 3E“3 , 281 . 22083 , 

1. 7191 75 ,0.4528E-3,0. ,23. 452294 ,-0.01 301 25, -0.1 64E-5/, 

RPD/. 01 74532925/, TWP/6. 2831 853/ , APT/. 446229 73E-8/ , VUPD/ 1 07 . 08 5/ 
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DECK 

SUND 


SUBROUTINE SUN 


12.01 

SUN 

4.27* 

12.09 

100 

12.10 

12.13 

110 

12.14 
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VI. SUBROUTINE MOON 


The relationships on page 47 and 48 of Reference 2 are used to compute 
lunar ephemeris quantities and their time derivatives (referred to the Earth's 
equatorial plane). (See page 23 of Reference 1.) 
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DECK MOOND 

C OHHON / OR B / 9 S ( 7 ) ,00S (T> , OL ( 7) , D9L ( 7) , A , £ ,01 » El , W » ZW f EN , I BUG , l , M 
DIMENSION C ( 1 4) 

0 A TA ( C ( 1 ) , I = 1 , 14) / .012150668 , *094900489 , 270. 43416,481267.883 
-0.00113,334.32956,4069.03403,-0.010325,-. 12E- 4, 259. 183275 , 
-1934.14201 , .00208,13.064993,5.1453964/ , 

RPO/. 01 74532925/ ,TWP/ 6. 283 1853/, APT/ .446229 73E-8/ ,VUPD/ 107. 08 5/ 

3300 FORMAT (1H021X5MSOLAR 17X3HDOT 15 X5HLUNAR17X3HDOT/ 

9X1HM4E20.6/9X1ME4E20.8/9X1MN4E20.8/4X6HLAMBDA4E20.8/ 
5X5HOMECA4E20.8/1X9MCAP OME6 A 4E 20 . 8 / 9 X 1 Hi 4E20 . 8 ) 
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DECK MOOND 


SUBROUTINE MOON 


15.01 

MOON 

4.284 

15.04 

10 

13.05 

13.15 

30 

13.16 

13,21 


13.18 
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COLE SATO 

SUBROUTINE NOON 


DECK MOOND 


LUNAR QUANTITIES 
PARAMETERS (T) 


SUBROUTINE MOON 


T2 = T 

♦ T 

TT = T 

♦ T 

«L(1> - 

C ( 1 ) 

QL ( 2 ) □ 

C (2) 

81(3) = C (13) * 

RPD/VUPD 


BEGIN OO LOOP 
10 1 = 1 , 3 
; 

SET I = 1 




10 

4 

D«L (I! 

1 = 0 . 


ADO 1 TO I , 
.COMPARE TO 3 


END Of 00 LOOP 
VARIABLE I 


QL<7> = 

ACOS (CLl ) 

SLI = 

SIN (QL (7) ) 

FAC 

= SW/SLI 

SLW = 

SIP * FAC 

GL(6) = 

ASIN (SLW) 

CZ = 

(CSI-C1P*CLI) / (SIP* 
SLI) 

z = 

ACOS (CZ) 

«L<5> = 

Z ♦ G - WP 

QL(4) 
«L (5) 

= QL ( 6) ♦ 
♦ CR - G 

10 


DQL(T) = 

(DQS ( 7) * (CIP*SSI* 
$1P*CSI*CW> -DW*SIP* 
SSI*$W> /SLI 



1 1 

OQL (6) = SIP * 
(DW*CW-FAC*CLl*OQL ( 
7 ) ) / (SLI4C0S (QL ( 6 ) ) 


12 

OQL (5) = DO - DW ♦ 
(DQS (7) *CSI*SW + DW* 
SSI*CW-SSI*CLI*FAC 
♦OQL (7) ) / ( SL 1 *C Z ) 


13 

OQL (4) = 

OQL ( 5 ) ♦ 

OQL (6) ♦ 
OCR - OG 


14 

BEGIN 
50 I 

OO LOOP 
= 4,7 
; 

SET 

1 = 4 





















VH. SUBROUTINE SRP 


The amplitude and phase of each term in the equations from pages 44, 45 
and 46 of Reference 2 for Solar Radiation Pressure is computed by this sub- 
routine. (See page 21 of Reference 1.) 


44 



OeCK SRFO 


0 1 MCN $1 ON RU3) ,c (5) 

common/ors/«S( 7) iDasrri , «i cri v o«L(n f A f c,oz v cl i tf,cv 9 eN v iftu««L tit 
CONMON/DCR/D (0) 'WO'CWO 
CON MOM /RAD/P ( • , 

t«UlVALCNCC(Sl,QS(?)),(SLO,OQS(4)),(5L,aS(4)) 
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DECK SRPO 
SUBROUTINE SRP 


14.01 

SRP 

5.064 

14.05 


14.08 

14.07 

100 


14.13 

200 

14.14 

14.20 


14.22 

14.21 

300 


14.28 

400 

14.29 

14.34 


14.36 

14.35 

500 


15.03 

600 

15.04 

15.09 


15.13 

15.12 

700 


15.18 

800 

15.19 

15.23 


15.27 

15.26 

900 
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COLE SAtD 


DICK SRPD 


PACE 14 


,j 


SUBROUTINE BtP- 


**PERTURBATlONS BY SOLAR 
RAOIATION PRESSURE** 
PARAMETERS (F) 


/ ^ 7 

s.o* 


SUBROUT1 

1 NE SRP 


2 


Twr 

* 1 

- - | 

SOI 

x 

SIN (01 ) 

SSI 

x 

SIN(SI) 

coi 

x 

COS (01) 

csi 

- 

COS(SI) 

HOI = 


* (1.-C0I) 

HSI = 


* (1.-CSI) 

VOI = 

*5 

* (l.*COI> 

VSI s 


* (1.+CSI) 

z - 


- E ♦ E 

RZ 

x 

SQRT(Z) 

R (?) 

= - 

• i .5 * r * 
E 

R (1) 

r 

R ( 7) * A 

R (2) 

= 

VOI * HSI 

R (3) 

= 

VOI * VSI 

R (4) 

X 

HOI * VSI 

R (3) 

= 

HOI * HSI 

R (6) 2 

- .5 * SOI * 
SSI 

R (0) 

= A/E 
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COLl SATO 


DEC A SR PC 


M6C 14 

(CONT) 


SOlROUTlNt **P 
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COLC \k TO 


DEC* |MD 


FACE 15 
\ (CONTI 


ftOttOUTlNK IM 



It 



END OF DO LOOP 
VARIABLE I 











Vm. SUBROUTINE SLGP 


The Solar and Lunar Gravitational Perturbations subroutine, SLGP, com- 
putes both long-period and secular effects using equations 4 and 5 from page 9 
of Reference 2. The amplitude of each of the long period terms is also com- 
puted (See page 20 of Reference 1) . 
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DEC * SLGPD 


DIMENSION 2 (67) , AC (23 ,7) 

C 1 MEN Si ON XAP(2) ,COEFA (8,2) 

COMMON /ORB/ 0 ( 14 , 2) , XA ,XE , XI , EL ,OMA ,COM , XN 

COMMON /PERT/ ARG ( 1 2 5 , 2) , ARGDOT (12 5,2) ,DE(125,2),DI(125,2) 1 DL(125, 
2) , OW ( 1 25 , 2) , DOM ( X 25 , 2 > , PER (12 5,2) , SEC (5 ,2) 

COMMON / C OE F / CLAMP (123) , COMAP (125) ,CCOMP (125) , COMA (125) ,CCOM (125) 
COMMON /DER/ D ( 6) , OM AD T , COMO T , E LD , 1 SEC ( 2 ) , K $L ( 2 ) 

real lamp.lamfdt 

DATA XAP(l) , XAP (2) /23440, ,60.2601/ 
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COLC satc 


DEC* 5L0PC 


TABLE OF CONTENTS AND CROSS-REFERENCE LISTING 


PACE SOX 

LABEL 

REFERENCES 

16.01 

SLGF 

4.31* 5.02* 

16.04 

110 


16.07 

120 

16.03 

16.08 

121 


16.09 

122 

16.07 

16.11 


16.13 

16.12 

10 


17.06 

801 

17.07 

17.10 

802 

17.11 

17.16 

803 

17.17 

18.01 

804 

18.02 

18.11 


18.13 

18.12 

910 


18.18 


18.20 

18.19 

911 


18.25 

915 

f 8 . 2 6 

19.02 

916 

19.03 

19.08 

917 

19.09 

19 . 17 


19.19 

19 .18 

918 


19.22 

919 

19.23 

19.26 


19.29 

19.28 

901 


19.32 


20.02 

20.01 

903 


20.06 


20.10 

20.09 

904 


20.15 


20.18 

20.17 

905 


22.42 


23.03 

22.44 


23.01 

22.45 

701 


23.05 

700 


23.06 


23.12 

23.08 


23.10 

23.09 

703 
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deck slgfd 
SUBROUTINE SLGP 


23.12 

7 * 0 ? 


23.15 


23.24 

23.17 


23.22 

23.21 

710 


23.24 

711 
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suMoutinc si tr / o \ 


\.5*SIS> * 
1 . 3*SI PS) 


ACS - 
Cl * 

- 4. * SI * 
(1 . -1 . 54SIPS) 

DUMB 

♦ 

(1 . 

= 2. * XNP**2 
XMP * XA * 
♦1.5*XE**2) 

AC 3 

= DUMB * AC3 

ACS = 

* 

DUMB * 3 . / 32 . 
XA * AC 5 

SEC (1 

.INDEX) = 0.0 

SEC (2 

.INDEX) = 0.0 

SEC (3 
• AC 3 

.INDEX) = CA4 
♦ CAS * AC4 

SEC (4,1 NDEX ) = - 

CA2 * AC 5 - C A 1 * 
AC4 

SEC ( 5 i I NDEX ) = - 

1 i 

A3 * AC 5 | 


END OF DO LOOP 
VARIABLE I 



C1P 

= COS(XIP) 

CIS 

= C I *42 

C1PS 

= C I P **2 

SSI 2 = 

SI N ( X I / 2 - ) **2 

SST P2 = 

SlN(XlP/2.)**2 

C SI 2 = 

CO$(Xl/2.)**2 

CSIP2 = 

C0S(XiP/2.)**2 

SF I 2 

- SSI 2**2 

SFIP2 

- SSI F2**2 

CFI 2 

= C SI 2**2 

CFIP2 

= C S I P 2**2 

S2IP = 

SIN(2.*XIP) 

S2I = 

SIN (2.*XI ) 

C 2 1 = 

COS ( 2 . *X I ) 

SI 2 = 

SI N ( X I / 2 . ) 

C I 2 = 

COS ( X I / 2 . ) 

CCI2 

= C I 2**3 

SC I 2 

= SI 2**3 

C6IP2 

= CSI P2**3 


THEP = 3. * EP 

Z (1) = 2. 

Z (2) = - EP 
Z ( 3 ) = SEP 
Z (4) = 2. 

7(A) r - THFP 

Z<6) = THEP 
Z ( 7) = 2. 

Z ( 8) = THEP 
Z ( 9 ) = - THEP 
Z (10) = 2. 

Z < 1 1 ) = - EP 
Z ( 1 2 ) = SEP 
Z (13) = 2. 

Z (14) r - EP 
Z ( 1 5 ) = SEP 
Z (16) = 2. 

Z (17) = - EP 

Z ( 1 8 ) = SEP 
Z ( 1 9 ) = 2. 

Z (20) = - EP 

Z C 21 ) 5 SEP 
Z (22) = 2. 

Z (23) = - EP 
Z ( 24 ) = SEP 
Z ( 2 5 ) = 2. 









COLE SATD 


DECK SlGPO 


03 / 23/67 


FACE 17 



z<«e> = - ep 

2(27) = SEP 

2(2®) = 2. 
2(29) = THEP 
2(30) = - THEP 
Z (31) = 2. 

2(32) = - THEP 
Z ( 33 ) = THEP 
Z ( 34 ) = 2. 
2(33) = - THEP 
2 (36) = THEP 
2(37) = 2. 

Z (38) = - EP 
Z ( 39) = SEP 
2(40) = 2. 

Z (41 ) = - EP 

2(42) = SEP 

2(43) = 2. 
2(44) = - EP 

2(43) = SEP 
Z ( 46) = 2. 

Z (4 7) = THEP 
2(48) = - THEP 
Z (49) = 2. 

Z(30) = THEP 


AC ( 3 , 1 ) = 4. 4 CFI2 
* SIPS 

AC (4,1) = - 8. * 

SFIZ * CFIP2 

AC (5,1) =-8.4 
CFI2 4 SFIP2 

AC (6, 1) = 8. 4 SFI2 
4 SFIP2 

AC ( 7 , 1 ) =3.4 SIS 
4 SIPS 

AC ( 8 , 1 ) =-3.4 
SIS 4 SIPS 

AC (9 , 1 ) = 8. 4 SI 4 
CSI 2 4 SIP 4 CSIP2 

AC (10, 1) = - 4. 4 
SI 4 CSI 2 4 S2IP 

AC ( 1 1 i 1 ) = - 4. 4 
SFI 2 4 SIPS 

AC (12, 1) = - 4. 4 

SI 4 SS12 4 S2IP 

AC ( 1 3 , 1 ) = 8. 4 SI 
4 SSI 2 4 SIP 4 
CSIP2 


AC ( 1 4 , 1 ) = 8. 4 SI 
4 CSI 2 4 SIP 4 
SSIP2 

AC ( 1 5 1 1 ) = 8. 4 SI 
4 SSI 2 4 SIP 4 
SSIP2 



DUMB = 2. 4 DUMB 

AC (1 ,2) = 4. 4 CFI2 
4 CFIP2 


AC ( 3 , 2) = 2. 4 CFI 2 

4 SIPS 

AC (4,2) = 4. 4 SFIZ 

4 CFIP2 

AC ( 5 , 2 ) = - 4. 4 
CFI 2 4 SF1P2 

AC (8,2) = - 4. 4 

SFI 2 4 SFIP2 

AC (7,2) =0.0 

AC (8,2) = 0.0 

AC ( 9 , 2 ) = 2. • SI 4 
CSI 2 4 SIP 4 CSIP2 

AC (10,2) = - SI 4 

CS12 4 S2IP 

AC (11 ,2) = 2. 4 
SFI2 4 SIPS 

AC (12,2) = SI 4 

SSI 2 4 S21P 

AC ( 1 3 , 2 ) =-2.4 
SI 4 SSI2 4 SIP 4 
CSIP2 




801 

6 

AC (I ,1) 

s 

o 

o 



14 

AC (14,2) = 2. 4 SI 
* C SI 2 4 SIP 4 
SSIP2 

AC (15,2) = - 2. 4 

SI 4 SSI2 4 SIP 4 
SSIP2 

1 

1! 

BEGIN 
803 1 

SET 

DO LOOP 
= 1, 15 
I 

* 1 
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E 1 S 


> V w w ll< 


subroutine slot 


AC (1 , 2) = AC (1 , 2) ♦ 

DUMB 




HI &H 

3 

END OF 

DO LOOP 

; VARIABLE 1 [ 


4 

1 AC (1,3) = 

.5 4 CF I 2 1 

4 CF1P2 

AC ( 2 , 3 ) = 

.25 4 SIS 

4 (1 . -1 

.54SIPS) 

AC (3,3) 

= .25 4 

CFI 2 

4 SIPS 

AC (4,3) = 

.5 4 SFI 2 

1 4 CFIP2 

I AC (5,3) = 

.5 4 CFI 2 I 

4 SFIP2 

j AC ( 6 , 3 ) = 

.5 4 SFI 2 I 

4 SFIP2 

AC (7,3) 

= 3./18. 4 

SIS 

4 SIPS 

AC (8,3) 

= 3 . / 1 6. 4 

SIS 

4 SIPS 

ACCt,3) = 

.5 4 Si 4 

C SI 2 4 SIP 4 CSIP2 

AC (10,3) 

= - . 25 4 j 

$1 4 C SI 2 4 S2IP 

AC (1 1 , 5) : <25 4 

1 * rtt 

4 SIPS j 
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DUMB = DUMB1 * 

3 . / 3 2 . * X A ! 

I 

DUMB1 = X5./64. * j 

XNF442 * XMF * i 
X A 442 * XE 442 j 

I 

AC ( 1 , 5) = - 8 • * 

SI 2 4 CC1 2 * CFI F2 i 

AC (2 , 5) = 4. * SI * 

Cl * (1.-1.54SIFS) 

AC (3,5) = -4.4 

S 1 2 * CC I 2 t SIPS 

AC (4,5) = 8. * C I 2 

* SC I 2 4 CFIP2 

AC (5 , 5) = - 8. * 

Si 2 * C C I 2 * SFIP2 

AC (6 , 5) = 8. * CI2 

* SC 1 2 * SFIP2 

AC ( 7 , 5) = 3. 4 SI * 

Cl * SIPS 

AC (8,5) = 3. * SI * 

Cl * SIPS 


AC (9 , 5) = 2. * 

( C 1 +C 2 1 ) * SIP * 
CSIP2 

AC (10, 5) - - 

( C l >C 2 1 ) * S21F 

AC <11 ,5) =4. * C12 

* SC 1 2 4 SIPS 

AC ( 1 2 , 5 ) = (CI-C2I) 
♦ S2IP 

AC (13 , 5) s - 2. * 

(C1-C21) * SIP * 
CS1P2 

AC (14 , 5) = - 2 % 4 
(CI+C2I) * SIP * 
SSIP2 

AC (15,5) = 2 . * 
(CI-C2I) 4 SIP 4 
SSIP2 



I 















COll S4TB 


DICK HMD 


M6C II 


SUBIOOTINC SLM 
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COlt SATO 


DECK 51 5PD 


PACK 1* 


SUBPOUT 1 HE SLOT 




11 


AC ( 1 , 7) = 4.4 CPI 2 

* CFIP2 

AC ( 2 i 7 ) = 0.0 

AC (3,7) = 2. * CFI 2 
* SIPS 

AC (4,7) = 4. * SF12 

* CFIP2 

AC { 5 , 7) = - 4. 4 

CFI 2 * SFIP2 

AC (6 , 7) = - 4. 4 

SFl 2 * SFIP2 

AC ( 7 , 7) = 0.0 

AC lb , n = u. u 

AC (9,7) = 2. 4 SI * 
C SI 2 4 SIP * CSIP2 

AC (10 , 7) = - SI 4 
C$12 4 S2IP 

AC ( 1 1 , 7 ) = 2. 4 
SFI 2 4 SIPS 

AC C 1 2 , 7) = SI 4 

SSI 2 4 S2IP 

AC (13,7) =-2.4 

SI * SSI2 4 SIP 4 
CSIP2 


12 


AC (14,7) = 2. 4 SI 
4 C S I 2 4 SIP 4 
SSIP2 

AC (15,7) = - 2. 4 
SI * SSI2 4 SIP 4 
SSIP2 

AC <16, 7) = - S2I 4 
S2IP 

AC (17,7) =-2.4 

SIS * SIPS 

AC (18,7) = - 4. 4 

SIS • CFIP2 


AC (19,7) = 2. 4 $21 

♦ 51 P * CSIP2 


AC (20,7) = 4. 

* SIS 

* SFIP2 


AC (21,7) * 2. 

* 521 

4 SIP 4 5S1P2 | 

AC (72 , 7) = 

0.0 

AC ( 2 3 , 7 ) = 

0.0 



0UM6 = 15./128. 4 
XNPM2/XN 4 XMP 4 
XE * SRE 


4 


(I -1 ) /3 ♦ 1 








ADO 1 TO l, 
COMPARE TO 23 


EMC OF DO LOOP 
VARIABLE I 





















tout sfcTo 


DECK SL&FC 


QJ/23/6 T 


FA&E 20 


SUBROUTINE SlfrT 
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coil UTD 


OtCK SLCPO 


P**E *0 



JCONT) 

0 


31 


de ( ?9, index) = 

s./s. * si * s t* 

SIP * $2 P 



32 

DE ( #0, INDEX) = - 

15. /4. * SIS * C SI 2 
* SIPS * CSIP2 


33 

DE ( SI, INDEX) * - 

15. /4. ® SI S • SSI2 
• SIPS • SSIP2 


34 

DE ( 62, INDEX) = 

5. /4. • $1 • SOR 4 
SCIP 


35 

DE( 83, INDEX) = - 

5 . / 4 . * C SI 2 * S 1 OR 
* SIPS * SSIP2 


36 

DE ( 84, INDEX) = - 

5 . / 4 . * SSI2 * 

SN10R * SIPS * 
CSIP2 


37 

DE ( 85,1 NDEX ) = - 

5./2. * SI * SN2R * 

sip * sr i p 2 


38 

DE ( 86, If 
5 . / 2 . * 
SIF * 

JDEX) = - 

>1 4 S2R * 
CFIP2 
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COLC n TO 


DECK SLCPD 


PA6C 


21 


subroutine SUSP 





67 
















coit mo 


mt ti 

ICQNT) 






























coie sue 


DEC*'. SLCFD 


PAGE 2 2 


subroutine slot 





69 









cod UTD 




Mk&C l l 
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COIC UTO 


DECK IKPO 


!1 

(CONT) 












deck bcdso 


DIMENSION A (Ml , MW) 
DIMENSION TEMP (20) 


DEC* BCD SO 


SUBROUTINE BCOSO 


24-01 

BCD SO 

6.05* 

24.04 

5 


24.05 

10 

24.33 

24.07 

15 


24. OB 


24.32 

24 .10 


24.14 

24.11 

20 


24.12 

77 

24.03 

24.14 

25 

24.10 

24.17 


24.18 

24.18 

40 


24.22 


24.26 

24.23 


24.24 

24.24 

45 


24.26 

50 


24.29 


24.30 

24.30 

55 


24.32 

60 


24.34 

75 

24.06 


7.24* 


75 









ADD 1 TO I Z , 
.COMPARE TO MW 


END OF DO LOOF 
VARIABLE 14 


LO/E y 

— <" ADD 1 TO IZZ, 

\comfare to m 

x / 


END OF DO LOOP 
VARIABLE I Z 


END OF DO LOOF 
VARIABLE IZZ 


















con i»n 

/^Saumoutinc *coso 


m 





nn 
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X. SUBROUTINE DATA 


This subroutine contains data only— this data is used by subroutines SRP , 
SLGP and the main program. The data in this subroutine is presented here. 
(See page 24 of Reference 1.) 
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DEC* OATAD 


COMMON/RAD/P (6, 6) ,B(6,6) 

COMMON /COEF / CLAMP(125) ,COMAP(125) ,CCOMP(125) ,COMA(125) ,CCOM(125> 
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COLE SATD 


DECK D A T A D 


FAC l 1 5 


THIS IS A BLOCK DATA DECK 

£Z7 



BLOCK DATA 


DATA (CLAMP (I ) , I = 

1 ,67)/2. 

,1- 

,3. 

,o. , 

rl. ,1. 

pO. 

pi . 

.l.i 

p2.p 

I.p3.p2.pl*p3.p2., 

p 1 . p 3 . p 2 . 

,1. 

,3. 

p2. , 

I 1 • 

p3. 

p2. 

pi. 

» 3. , 

,0. , 

l.pO.pl.pl. pO.pl*: 

1 1 . , 2 . p 1 . 

,3. 

p2. 

. 1 . . 

,3. 

,2. 

.»• 

»3. 

pO. 

,1. p 

1. ,1. ,2. ,1. ,3. ,2. 

p 1 • i 3 • p 2 . 

,1. 

,3. 

,2- 

pi. 

,3. 

,2. 

pi. 

,3. 

pi./ 

DATA (COMAP ( 1 ) , 1 = 

1 , 67) /O. 

,1. 


. ,0, 


-1 . 

,0. 

p-1 

‘.-1 


0. | 1. ,-1 . |0. , 1 . f -l . | 0 • 1 1 • p ** 1 • p 0 . pi. p * 1 • p 0 . pi. , - 1 . , 0 . pi. i*l . p 
0 • t * 1 * i '1 ■ »0. i * 1 • p * 1 * pO« i *1 * t “1 • P o • i 1 ■ > “1 • »0. 1 1 • i • l • ,0. p 

1 . , - 1 . , o. , - 1 . ,-l. , 0 . , - 1 . , - 1 . , 0 . pi. ,-l. , 0 . pi. ,- 1 ., 0 ., 1 .,-t., 

0 . .1. f -l., 0 .,l.,*l.,-l./ 

DATA (CCOMP ( I ) , I = 1 ,« 7 )/ 0 ., 1 . ,-l., 0 .,-l. ,-l. , 2 . , 1 . ,- 3 . , 0 . , 

1 . p*l. p* 4 . p- 3 . p* 3 « p* 4 * p - 3 . »* 3 . i - 2 . p * 1 • >* 3 . >* 2 . p*l * > * 3 ■ » 

* 1 . i 0 . p - 2 . pi. p 0 . p * 2 . p 2 . p- 3 . pi. pi. p - 2 . pO. p*l. pO. p * 2 . p* 3 . t 

- 2 . f - 4 . ( - 3 .,- 2 .,- 4 . pl.pO. p- 2 .p 2 .pl. ,* 3 . p 0 . pi. ,-l.p-l. t 0 . p- 2 . , 

* 4 . p * 3 . p - 5 . p- 3 . f - 2 . p - 4 . p * 2 . p “ 1 • p “ 3 . p*l./ 

DATA (COMA (I) , I = 1 , « 7 ) / * 2 . , * 2 . , - 2 . , 2 . , - 2 . , 2 . , - 2 . , - 2 . , 

2. p 2 . p 2 • i 2 ♦ p 2 . p 2 . p 2 • p*2. p~2. p * 2 . p*2. p*2. p*2. t 2. p 2 . , 

2. p*2. p*2. p-2. , -2 . p*2. p 2 . f 2. p*2. p2. , 2 • p*2. p 2 • $ 2 • p 2 . p 

2 . . 2 . . 2 . . 2 . , - 2 . p- 2 . p- 2 ., 0 ., 0 ., 0 ., 0 ., 0 ., 0 ., 0 ., 0 ., 0 ., 0 ., 0 ., 0 ., 0 ., 

O.pO.pO.pO.pO.pO.pO.pO.pO. / 

DATA (CCOM ( 1 ) , I = 1 , « 7 ) /* 2 . , - 2 . , - 2 . , 0 . , 0 . , 0 . , - 2 . , -2 . , 2 . , 

- 2 . f - 2 .,- 2 . / 2 ., 2 ., 2 ., 2 .p 2 ., 2 ., 0 ., 0 ., 0 ., 0 ., 0 .p 0 .p-l.p-l.p-l., 

- 1 . p- 1 . pi. p* 2 . p 2 . p * 2 . p-l.pl.p-l.p-l.p-l.p- 1 . pi. pi. pl» pi. pi. p 

DATA ( (B ( 1 > J > , 1 * 1 , 6 ) ,J«I, »“*• 
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COLC SAID 


OAT.A (CCQN(l) ,1 = 1 ,6?>/-2. ,-2. ,-2. ,0. ,0. ,0. ,-2. ,-2. ,2. , 
-t.»-t.,-t.it. > t.,t.i't.|t.,t. l 0.|0.>0.|0..0. l 0.|-t..-t..-l.i 
-1. .-1. ,1. ,-t. ,2. ,-2. ,-l. • I. ,-l. ,-1. ,-l. ,-l. ,1. .I. it. ,1. ,1 
l.»-l..-l.it.i-t.,-t.|t.,-t.,-t.,-t.|-l..-l.,-l.,t.it.,t'it.,t 

1 . ,0. ,0. ,0. ,0./ 

DATA (CLAMP tl) ,1=68,123) / 14*1 . 0 , 1 4*3 . 0 , 1 4*1 .0,1463.0,280.0/ 

DATA (COMAP (1 ) , l = 68 , 1 2 5 ) / 5S*0 . 0 / 

DATA (CCOMP(l) ,1=68 ,125) /-l. 0,-2. 0,0. 0,-3. 0,1.0, -4.0, 

2. 0, -1.0, 0.0, -2. 0,1. 0,-3. 0,2. 0,-4. 0,-3. 0,-4. 0,-2. 0, 

-5. 0,-1. 0,-6. 0,0. 0,-3. 0,-2. 0,-4. 0,-1. 0,-3. 0,0.0, -6.0, 

-1.0, -2. 0,0. 0,-3. 0,1. 0,-4. 0,2. 0,-1. 0,0. 0,-2. 0,1.0, -3.0, 

2. 0, -4. 0,-3. 0,-4. 0,-2. 0,-5. 0,-1. 0,-6. 0,0.0, -3.0, 

-2. 0,-4. 0,-1 . 0,-5. 0,0. 0,-6. 0,0. 0,0.0/ 

DATA (COMA (1) , 1 = 68 , 1 25 ) / 7*1 . 0 , 7*- 1 . 0 » 7*1 . 0 , v 
7*-l .0,7*3. 0, 7*-3. 0,7*3. 0,7*-3. 0,0. 0,0.0/ 

DATA (CCOM ( I ) , 1=68, 125 )/0. 0,1. 0,-1. 0,2. 0,-2. 0,3.0, 

-3. 0,0. 0,-1. 0,1. 0,-2. 0,2. 0,-3. 0,3. 0,0. 0,1. 0,-1. 0,2.0, 

-2. 0,3. 0,-3. 0,0. 0,-1. 0,1. 0,-2. 0,2. 0,-3. 0,3. 0,0.0, 

1.0, -1.0, 2. 0,-2. 0,3. 0,-3. 0,0. 0,-1. 0,1. 0,-2. 0,2.0, 

-3. 0,3. 0,0. 0,1. 0,-1. 0,2, 0,-2. 0,3. 0,-3. 0,0. 0,-1. 0,1.0, 

-2. 0,2. 0,-3. 0,3. 0,0. 0,0.0/ 
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ENTRY 


ROBIN 



DECK DATAD 


DECK RDB 


26.01 

ROBIN 

1.22* 

26.02 

ST 

26.02 

26.04 


26.03 

26.05 


26.04 

26.07 


26.07 

26.12 

ROC 

26.08 

26.15 

ERR 

26.09 

26.17 

RNC 



26.14 


26.13 
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